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Abstract 

We construct a supersymmeterized version of the model presented by Grimus and Lavoura 
(GL) in [1] which predicts O23 maximal and #13 = in the lepton sector. For this purpose, we 
extend the flavor group, which is D± x z^ 1 in the original model, to D4 x Z5. An additional 
difference is the absence of right-handed neutrinos. Despite these changes the model is the 
same as the GL model, since O23 maximal and 6*13 = arise through the same mismatch of D4 
subgroups, D2 in the charged lepton and Z2 in the neutrino sector. In our setup D4 is solely 
broken by gauge singlets, the flavons. We show that their vacuum structure, which leads to the 
prediction of 13 and #23, is a natural result of the scalar potential. We find that the neutrino 
mass matrix only allows for inverted hierarchy, if we assume a certain form of spontaneous CP 
violation. The quantity |m ee |, measured in neutrinoless double beta decay, is nearly equal to the 
lightest neutrino mass m^. The Majorana phases <f>\ and 4>2 are restricted to a certain range for 
m 3 ^5 0.06 eV. We discuss the next-to-leading order corrections which give rise to shifts in the 
vacuum expectation values of the flavons. These induce deviations from maximal atmospheric 
mixing and vanishing #13. It turns out that these deviations are smaller for 623 than for #13. 
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1 Introduction 



Experiments revealed that neutrinos have properties very distinct from the other known fermions. 
Their masses are much smaller and their hierarchy is much less pronounced. Unfortunately, only 
two mass squared differences, the solar and the atmospheric one, are known from experiments [2] 

Am 2 21 = (7.65±g;|g) • HT 5 eV 2 and \Am 2 31 \ = (2.40±g;||) • 10~ 3 eV 2 (2a) (1) 

where Am^ denotes mf — m 2 with rriij being the neutrino masses. As the sign of Am^ is unknown 
it is not clear whether neutrinos follow a normal or an inverted hierarchy. Even more surprising 
than the neutrino masses is the fact that leptons have a very peculiar mixing pattern with two 
large mixing angles and one small one. The experimental measurements [2] 

sin 2 6>i3 < 0.040 , 6»i3 < 11.5° , 0i 3 < 0.2 , 

sin 2 23 = O.5O^? , 23 = (45.0+™)° > 023 = 0.7851°;^ , (2a) (2) 

sin 2 0x2 = 0.304^6 ? 0i2 = (33.5+2.8)0 ? 12 = 0.58^ , 

allow for the possibility that the atmospheric mixing angle 023 is maximal and the reactor mixing 
angle 0i3 vanishes. These values can be deduced from a neutrino mass matrix which is /i — r 
symmetric [3], i.e. does not change its form if second and third columns and rows are interchanged, 
in the charged lepton mass basis. At the same time the solar mixing angle 0i2 is left undetermined. 
An even more constraining pattern is the one of tri-bimaximal (TB) mixing [4] in which, apart 
from 023 = 7r/4 and 0i3 = 0, also 6 12 is fixed by sin 2 0i2 = 1/3. Both patterns have been subject to 
extensive studies in order to find a theoretical explanation. The most promising one seems to be 
the assumption of an additional flavor symmetry which is responsible for such a mixing pattern. 
For TB mixing the simplest models are based on the tetrahedral group A4 [5], while for n~ t 
symmetry even smaller groups are appropriate such as the dihedral groups D3 = S3 [7] and D4 [1]. 

An interesting observation which has been made first in the models and then also in the 
models predicting /i — r symmetry with a dihedral flavor group is the fact that the VEVs of a certain 
subset of scalar fields {4> e }, coupling only to charged leptons at leading order (LO), preserve one 
subgroup of the original symmetry, while another set of scalars {<j> v }, coupling only to neutrinos at 
LO, breaks the flavor group to a different subgroup. This mismatch can be regarded as an intuitive 
explanation for sizable mixings in the lepton sector. For more general considerations on the 
origin of a certain mixing pattern see [8]. 

One of the main issues in these models is the vacuum alignment of the flavor symmetry breaking 
fields, since without a special alignment the mixing pattern is merely a result of parameter tuning. 
In a class of models [12] this alignment is achieved in a supersymmetric framework where the 
scalars transforming non-trivially under flavor are only gauge singlets. Their vacuum expectation 
values (VEVs) are driven by another set of gauge singlets, the driving fields. This mechanism to 
align the flavon VEVs is also used in models with other discrete [13] and with continuous flavor 
symmetries [14]. 

The D4 model [1] by Grimus and Lavoura (GL), which successfully predicts fi — r symmetry in a 
natural way including a profound explanation for the vacuum alignment, is non-supersymmetric in 

From a group theoretical point of view S4, the permutation group of four distinct objects, might be even more 
appropriate [6]. 

2 A similar observation has been made concerning the Cabibbo angle in the quark sector whose value might be the 
result of a non-trivial breaking of a dihedral group such as D7 or D14 [8-10]. 

3 Notice that there are also models which can predict lepton mixings without preserving non-trivial subgroups in 
all sectors of the theory, for instance [11]. 
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its original version. For several reasons, it would be desirable to supersymmeterize this model. In 
doing this, it is of advantage to break the flavor and the gauge group (spontaneously) by separate 
sets of scalar fields, flavons and Higgs doublets. In this paper we present such a supersymmeterized 
version in which the vacuum alignment is achieved in a similar way as in [12]. We arrive at 
#23 = 7r/4 and #13 = at LO and analyze the next-to-leading order (NLO) effects which will 
perturb this result in a particular way. In the minimal supersymmetric extension presented here 
the flavor symmetry D4 is accompanied by a Z5 symmetry which plays a similar role as z^ ux ^ in 
the original model. Furthermore, our model does not incorporate right-handed neutrinos so that 
the light neutrino masses stem from the dimension-5 operator lh u lh u /A. Despite these changes the 
model is essentially a supersymmeterized version of the GL model, since the prediction of maximal 
623 and vanishing #13 is still due to the fact that we preserve, at LO, a D2 subgroup of D4 in the 
charged lepton and a Z2 subgroup in the neutrino sector. Since this Z2 group is not contained in 
the D2 group of the charged lepton sector, D4 is completely broken in the whole theory. Apart 
from predicting the value of #13 and 623, the original GL model also predicts that neutrinos are 
normally ordered and that the effective Majorana mass of neutrinoless (3(5 (Qv(3/3) decay is equal 
to \m ee \ = 77117712/7713. These predictions result from the fact that in the GL model not all fields 
are present which are allowed to have a non-vanishing VEV in accordance with a preserved Z2 
subgroup in the neutrino sector. Here we will include all flavons allowed by the symmetry principle 
so that we still predict 623 = vr/4 and #13 = 0, but now can accommodate both mass hierarchies. 
In this respect our results are analogous to another model by GL predicting #13 and 623 with the 
help of the dihedral group D3 = S3 [7] . 

In our phenomenological study we concentrate on the possibility of a certain type of spontaneous 
CP violation in order to make the model more predictive. We find that neutrinos then have to have 
an inverted hierarchy and \m ee \ ~ 7773 holds. Furthermore, the Majorana phases 4>i t 2 can only take 
values in a limited range for 7773 < 0.06 eV. If we additionally remove one of the flavons from the 
model (this is analogous to what is done in [1]), the three parameters of the model in the neutrino 
sector are determined by the three measured quantities, Amfi, (Am^l and 612, and the Majorana 
phases are predicted to be <j)\ = ir/2 and 02 = 0. The NLO corrections coming from the inclusion 
of operators with one more flavon lead to deviations from 623 = vr/4 and #13 = 0. It turns out that 
#23 — 7r/4 is much smaller than #13. 

In [15] it has already been attempted to build a -D4 model in which the Higgs doublets trans- 
forming non-trivially under the flavor group are replaced by flavons. Since this model is non- 
supersymmetric the vacuum alignment problem is not straightforward to solve and indeed one has 
to require that one of the quartic couplings in the potential vanishes. However, such an assumption 
will not be stable against corrections and has to be considered as a severe tuning. In a second 
version of this model [16] which is supersymmetric, the potential is not studied such that the ques- 
tion of the vacuum alignment also remains open. Hence, a successful supersymmeterization of the 
original D4 model by GL still does not exist. 

The paper is organized as follows: in Section [2] we repeat the necessary group theory of -D4 and 
the properties of the subgroups relevant in the D4 model. Section [3] contains the LO results for 
the lepton masses and mixings as well as the flavon potential. In the following section the NLO 
corrections are studied. We summarize our results and give a short outlook in Section [5j In the 
two Appendices we treat additional group theoretical aspects of the model. 
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Table 1: Character table of the group D4. d are the classes of the group, °d 
is the order of the i th class, i.e. the number of distinct elements contained in this 
class, °hc i is the order of the elements S in the class &, i.e. the smallest integer 
(> 0) for which the equation S c * = 1 holds. Furthermore the table contains one 
representative for each class d given as product of the generators A and B of the 
group. 



2 Group Theory of D4 

In this section we briefly review basic features of the dihedral group D±. Its order is eight, and it 
has five irreducible representations which we denote as lj, i = 1, ...,4 and 2. All of them are real 
and only 2 is faithful. The group is generated by the two generators A and B which can be chosen 
as [17] 

a-(; °) ™.b-(;j) ( 3, 

for 2. Note that A is a complex matrix, although 2 is a real representation. For {a\,a^) T ~ 2 
therefore (a^, cl\) T transforms as 2 under D±. The generators of the one-dimensional representations 
can be found in the character table, displayed in Table [TJ The generators fulfill the relations 

A 4 = 1 , B 2 = 1 and ABA = B . (4) 

The product rules for 1: are the following 

li x li = ll . lj. X lj = for i = 1, 4 , 1 2 X I3 = I4 , 1 2 x 1 4 = 1 3 and 1 3 x 1 4 = 1 2 . 

For Si ~ lj and {a\,<Pi) T ~ 2 we find 

( Siai ) ~ 2 , ( S2ai ) ~ 2 , ( S3fl2 U 2 and ( ^ )~2. 

The product 2x2 decomposes into the four singlets which read for (oj, a,2) T , (b±, 62 ) T ~ 2 
ai&2 + 0-2^1 ~ ii j ai^2 — «2&i ~ i2 1 + °2^2 ~ 1^3 arid ai&i — 0262 ~ I4 • 

More general formulae for generators, Kronecker products and Clebsch Gordan coefficients can be 
found, for example, in [9,18]. Notice that our group basis does not coincide with the one chosen 
by GL in [1]. Therefore, the mass matrices shown below have another appearance, especially 
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the charged lepton mass matrix is not diagonal in our basis. However, the prediction of the 
mixing angles does not depend on the chosen group basis. In Appendix |A] we explicitly discuss the 
correlation between our basis and the one found in [1]. 

All subgroups of are abelian: Z2 = D\, Z± and D2 = Z2 x Z2. We are interested here in Z2 
subgroups which are generated by B A m with m = 0, ...,3 and the D2 subgroup generated by A 2 
and BA. In order to see that B A m gives a Z2 group note that 

(BA m ) 2 = BA m B A m = BA m_1 B A m_1 = • • • = B 2 = 1 

holds, if Eq.(|l]) is used. Similarly, one finds for A 2 and BA 

(A 2 ) 2 = A 4 = 1 and (BA) 2 = BABA = B 2 = 1 

by using again the generator relations. Obviously, A 2 and BA are not equal (in general) and thus 
they generate different Z2 subgroups. Additionally, we have to check that A 2 and BA commute 

A 2 BA = A 3 BA 2 = A 4 BA 3 = BAA 2 . 

All this shows that A 2 and BA generate a Z2 x Z2 group which is isomorphic to a D2 group. The 
other non-trivial element of the D2 group is BA 3 . Thus, one could also use the two elements A 2 
and BA 3 to generate this D2 group. However, we follow the convention to use as generators A 2 
and the element BA P with p being the smallest possible natural number. The Z2 symmetry given 
through BA m is left unbroken by a non- vanishing VEV of a singlet transforming as 1% if m is even 
and of one transforming as I4 for m being odd. Additionally, it is left intact by fields ^1,2 forming 
a doublet, if their VEVs have the following structure 

( IS )«('?)■ 

For preserving the D2 group generated by A 2 and BA only singlets in 1a are allowed to have a non- 
vanishing VEV. Especially, no fields forming a doublet under D4 should acquire a VEV. Clearly, in 
all cases singlets in the trivial representation of D4, lj, are allowed to have a non- vanishing VEV. 
Note also that in none of the cases a field transforming as 1.2 can acquire a non-zero VEV. Since we 
concentrate on the D2 subgroup induced by A 2 and BA, the Z2 subgroup has to be generated by 
BA m with m being even in order not to be a subgroup of the D2 group. Only then the mismatch 
between the two subgroups is achieved. The choice of m, m = or m = 2, depends on the relative 
sign between (^i) and (^2) for two fields -01,2 ~ 2. 

3 The Model at Leading Order 

We augment the Minimal Supersymmetric Standard Model (MSSM) by the flavor symmetry D4 x 
Z5. As mentioned above, the non-trivial breaking of D4 is responsible for maximal atmospheric 
mixing and vanishing #13, while Z^ is necessary to separate the charged lepton and the neutrino 
sector. The model contains three left-handed lepton doublets 1%, the three right-handed charged 
leptons e£, the MSSM Higgs doublets h u> d and two sets of flavons {x e , Ve} and {xv, <Pv, ^1,2} which 
break D4 in the charged lepton and the neutrino sector, respectively. The transformation properties 
of these fields are collected in Table El 
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Table 2: Particle content of the model. 
SU(2)l doublets, ef are the right-handed 
Higgs doublets. The flavons Xe, ty?e, Xv, fi 

2-Ki 

The phase factor wise s . 



h denotes the three left-handed lepton 
charged leptons and h Ut d are the MSSM 
and lpi $ 2 only transform under D4 x Z5. 



3.1 Fermion Masses 

The invariance of the charged lepton and neutrino mass terms under the flavor group D4 x Z5 
requires the presence of at least one flavon. Thus, charged lepton masses are generated by non- 
renormalizable operators only. In a model which treats quarks as well this allows the explanation 
of the small r mass compared to the top quark mass without relying on a large value of tan (3 = 
{h u )/{hd) = Vu/vd- The neutrinos receive Majorana masses through the dimension-5 operator 
lh u lh u /A which can be made invariant under the flavor group by coupling to a flavon. The part of 
the superpotential giving lepton masses reads at LO 

•W = ylXehe\^+ylXe(hel + hel)^+yy e (l 2 el-he^ (6) 
h 2 h 2 h 2 h 2 

h 2 

+ViXu{hh + hh)-j^- 

A is the cutoff scale of the theory whose order of magnitude is determined by the scale of the light 
neutrino masses, see below. For the moment we assume that the flavons Xe and (p e acquire the 
VEVs 

(ip e ) = u e and (xe) = w e . (7) 

As discussed in Section [2] these VEVs break D4 down to D2 generated by A 2 and BA in the charged 
lepton sector. The VEVs of the flavons coupling only to neutrinos at LO, are of the form 

(<Pu) = u , ( X u) = ID , ( ) = v ( \ ) , (8) 

and therefore leave a Z2 subgroup, generated by B, unbroken. As mentioned, the equality of the 
VEVs of {ijji) and {1P2) is crucial. As will be discussed in Section 13.31 the vacuum structure in 
Eq.© and Eq.([8|) is a natural result of the minimization of the flavon potential. We obtain the 
following fermion mass matrices, when inserting the flavon VEVs and {h u ^) = v u ^ 

Vd ( V ^ We ° ° \ v 2 ( VlW V2V V2V \ 

M l = ~T \ 2/3 u e V2 w e and M v = I y 2 v y 3 u y A w . (9) 

\ y e 2 w e -y%u e ) \ y 2 v y A w y 3 u J 

Thereby, the left-handed fields are on the left-hand and the right-handed fields on the right-hand 
side for M\. The matrix M\ Mj is diagonalized through the unitary matrix Ui, i.e. uj Mi M} Ui is 



5 



diagonal. \J\ acts on the left-handed charged lepton fields and is given by 




Ui = e i7r / 4 /V2 e- i7r / 4 /V2 . (10) 



For the masses of the charged leptons we find 

"i-e = -jr|2/iU>e| = \y 3 u e + iy 2 w e \ and m r = — \y 3 u e - iy 2 w e \ . (11) 

In order to arrive at non-degenerate masses for the \i and the r lepton either y%u e or y\w e has to 
be non-real indicating CP violation in the Yukawa couplings and/or flavon VEVs. For m T being 
around 2 GeV we find that for small tan (3 - corresponding to of the order of 100 GeV - the ratio 
of the flavon VEVs u e and w e over the cutoff scale A should fulfill @ 

^,^A 2 «0.04 (12) 
A A 

with A being the Cabibbo angle. The smallness of the ratio m e /m T is in this model only explained 
by the assumption of a small enough coupling yf. Similarly, m^/m T enforces a certain cancellation 
between the two contributions y%u e and iy\w e in m^. In [1] these problems have been solved by 
the assumption that the electron couples to a Higgs field different from those coupling to the [i and 
the r lepton and by an additional symmetry which leads to = 0, if it is unbroken, respectively. 
The neutrino mass matrix in the charged lepton mass basis reads (indicated by a prime ( ' )) 



2 / Viw y2V y2V 
A 2 



M' v = UjM v U[- = ^|| y 2 v Vi w y 3 u | . (13) 

y 2 v y 3 u y±w 



As M' u is /i — r symmetric, it immediately follows that the lepton mixing angle #13 vanishes and 
#23 is maximal. The solar mixing angle 6\ 2 is not predicted, but in general expected to be large. 
Also the Majorana phases 4>i,2 are not constrained. The lepton mixing matrix is of the form 



Umns = diag(e 



»7l 



/ cos 9 12 

-^T T "72 I ■diag(e^V^ ) . (14) 

sin 0i2 
V V2 



sin 9 12 





COS 012 


1 


V2 


V2 


COS 012 


1 


V2 


v/2 



The Majorana phases <f>i^ = 01,2/2 can be extracted from Umns by bringing it into the standard 
form [19]. Due to the additional factor 1/2 the phases </>i2 vary between and ir. Assuming that 
all flavon VEVs are of the same size, the estimate in Eq. f)12j) also holds for the VEVs of the flavons 
Xv, and ^1,2- For small tan/3, i.e. v u w w 100 GeV, a light neutrino mass scale between 
yj\ Am^l ~ 0.05 eV and 1 eV fixes the range of the cutoff scale A to be 

4 • 10 11 GeV < A < 8 • 10 12 GeV . (15) 

As shown in Section 13.31 we can assume that CP is only spontaneously violated in this model 
by imaginary VEVs w e and w of \e and %v- Thus, apart from w e and w all other parameters, i.e. 
couplings and VEVs, are real in the following. According to Eq. ljlip an imaginary w e allows the 



4 Although not excluded, there is no obvious reason to assume that there is a large hierarchy among the different 
flavon VEVs. In general, these are correlated through the parameters of the flavon potential. 
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\x and the r lepton mass to be non-degenerate. In the neutrino sector only the VEV w of \v is 
imaginary, whereas all other entries in M' v are real, so that the matrix in Eq. (]13p can be written as 



v 2 v 



i s 


t 


t 


t 


i x 


z 


t 


z 


% X 



where we define the real parameters 

Im(u>) Im(w) u . s 

s = yi , t = y 2 , x = y 4 and z = y 3 - . (17) 

V V V 

3.2 Phenomenology 

In the following we analyze the phenomenology of this model. For the eigenvalues of M'^M'J we 
find 



m 



2,1 



1 / vl 



2 / \ 2 



s 2 + At 2 + x 2 + z 2 ± ^{s - x) 2 (8t 2 + {s + x) 2 ) + 2(4t 2 + x 2 - s 2 )z 2 + z 4 



;2 



2 / \ 2 



and mi = (x* + z 2 ) . (18) 

This assignment of the eigenvalues is unambiguous, since m\ > m\ is experimentally known and 
the eigenvalue corresponding to the eigenvector (0, 1, — 1) T can only be m\. The solar mixing angle 
#12 is found to depend on s, t, x and z in the following way 

x z + z z — s z 

Before discussing the general case with unconstrained parameters s, t, x and z we comment on 
the special case in which z vanishes, since then the model contains three real parameters which can 
be determined by the three experimental quantities Am^, lAm^jJ and 9\2- According to Eq. (|17p 
either 1/3 or u have to vanish for z = to hold. Assuming that 7/3 is zero however has to be regarded 
as fine-tuning. In contrast to that, a vanishing VEV u can be explained either through the absence 
of the flavon ip v from the model or through a flavon potential which only allows configurations with 
u = to be minima. The neutrino mass is then proportional to \x\. From Eq. (|18p and Eq. (|19p 
we can derive for z = 

m2 = _ 1 cos 4 9 12 (Am 2 21 + Am| 1 (tan 4 j7i2 - l)) 2 

1713 4 sin 2 6»i2 Amjj^l + tan 2 #i 2 ) - Am 2 x ' ^ ' 

Neglecting the solar mass squared difference we can simplify this expression to 

ml ps -Am§! cot 2 26*12 • (21) 

Eq. (|2ip shows that Am^ < 0, i.e. the neutrinos have to have an inverted hierarchy. Note that 
similar results can also be found in [20]. A relation analogous to Eq. ([20|) can be found for \m ee \ 
measured in 0u(3(3 decay experiments. Note that \m ee \ is proportional to \s\ due to Eq. (fT6l) and 
can be written in terms of 7713, tan #12 and the mass squared differences as 

,2 _ m 2 - 2 tan 2 gig) + Amj^tan 4 9 12 - l)) 2 , , 

1 ee| " 3 (Am^ + Aml^tan 4 ^-!)) 2 ' 1 ' 
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Figure 1: \m ee \ plotted against for z ^ 0. The 
dashed (red) line indicates the results for z = 0. Mass 
squared differences and the solar mixing angle are in 
the allowed 2a ranges [2]. As one can see, \m ee \ and 
mz have nearly the same value. Additionally, one 
finds that has a lower bound around 0.015 eV. 
For z = we also find an upper bound on ms.. 



Figure 2: tan 6\2 plotted against mz for non- 
vanishing z. Again the dashed (red) line indicates 
z = (assuming the best-fit value for the atmospheric 
mass squared difference) and gives a lower bound for 
Apart from that the results for tan #12 are only 
constrained by the requirement that they are within 
the experimental 2<r ranges [2], 0.61 < tan 612 < 0.73. 



In the limit of vanishing solar mass splitting we find 

I "lee I ~ m 3 • (23) 

Taking the best-fit values Am| = 7.65 • ICT 5 eV 2 ,Am^ = -2.40 • 1CT 3 eV 2 and sin 2 12 = 0.304 [2] 
we obtain s « 0.02075, t « 0.03502, x 0.02146 | for v u 100 GeV, A « 4 ■ 10 11 GeV and v/Arn 
A 2 0.04. The neutrino masses are mi m 0.05348 eV, m 2 0.05419 eV and m 3 0.02146 eV. 
Their sum Yl m i ~ 0.1291 eV lies below the upper bound required from cosmological data [21]. 
|m ee | equals 0.02075 eV which might be detectable in the future [22]. The two Majorana phases 
01.2 are <j)\ = n/2 and $2 = 0. For tritium (5 decay we find mp ~ 0.05370 eV which is about a 
factor of six smaller than the expected sensitivity of the KATRIN experiment [23] . 

Turning to the general case with z/Owe first observe that also in this case the light neutrinos 
have to have an inverted hierarchy. To see this let us assume that the matrix in Eq. (|16p would 
allow the neutrinos to be normally ordered, i.e. 711,3 > m i as we h as m 3 > m 2- From m| — m 2 > 
then follows 

x 2 + z 2 - s 2 - At 2 - s/(s - x) 2 (8t 2 + (s + x) 2 ) + 2(4t 2 + x 2 - s 2 )z 2 + z 4 > . (24) 
From this we can deduce 

x 2 + z 2 >s 2 + 4t 2 and I6t 2 (t 2 + x(s - x) - z 2 ) > . (25) 
Rearranging the first inequality and taking t 7^ (otherwise #12 is zero) for the second one, we get 

x 2 - s 2 > At 2 - z 2 and t 2 - z 2 > x(x - s) . (26) 

The sum of these inequalities leads to 

5 Actually we find four solutions which all lead to the same absolute values, but to different signs for s, t and x, 
with the constraint that s and x have the same sign. 
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Figure 3: The Majorana phases (j>i (blue/darker 
gray) and <j>2 (green/lighter gray) plotted against the 
lightest neutrino mass mz for non- vanishing z. The 
values for z = 0, 4>\ = 5-, cj>2 = 0, are displayed 
by dashed (red) lines. Notice that the results for 
z / are centered around these values. The mea- 
sured quantities, Am^i, |Am§i| and 612, are within 
the 2cr ranges [2]. 



Figure 4: Phase difference cj>i — cf>2 against 7713 for 
z / 0. The case z = 0, |0i — <f>i\ = tt/2, is 
given by the dashed (red) lines. As one can see, 
0i — <j>2 1 is restricted to the interval [7r/2, 37r/4] for 
m 3 ^ 0-06 eV. Its deviation from 7r/2 increases with 
increasing 7713. Again, the mass squared differences 
and 612 are within the experimentally allowed 2<r 
ranges [2]. 



s(x - s) > 3t 2 > . (27) 

From Eq. (|27[) we see that s and x have the same sign, while x 2 > s 2 , hence x(x — s) > s(x — s). 
Combining Eq. (|26p and Eq. (|27|) . we find t 2 — z 2 > 3t 2 , an obvious contradiction. Thus, the neutrinos 
cannot be normally ordered as assumed by m 2 > m 2 . Instead we always have to 2 . > m| which is 
only possible in case of an inverted hierarchy. Note that it is a priori not clear that also mi is larger 
than TO.3, since the size of the mass squared differences has to be tuned so that Am 2 ^ <C lAmjjjJ. 
In fact, Am^i is given by 

= (x) 2 (J) 2 V(s - a0 2 (8t 2 + (s + x) 2 ) + 2(4* 2 + x 2 - s 2 )z 2 + z 4 . (28) 

It vanishes, if z = and s = x. Thus, Am^ <C (Am^l holds, if these equalities are nearly met. 
As noted, the vanishing of z can be made a natural result of the model. The near equality s ~ x 
however has to be regarded as a certain tuning of the couplings y% and j/4, see Eq. fjTT]) . 

We study the general case with a numerical analysis. To fix the light neutrino mass scale 

we adjust the resulting solar mass squared difference to its best-fit value. At the same time the 
atmospheric mass squared difference and the mixing angle #12 have to be within the allowed 2a 
ranges [2]. First, we note that our numerical results confirm that z has to be in general smaller 
than the parameters s, t and x and that s and x have to have nearly the same value. In Figure [T] we 
plotted \m ee \ against the lightest neutrino mass rn.3. As one can see, the approximate equality of 
\m ee \ and TO3, deduced for z = in Eq. (|23p . still holds for z 7^ 0. The dashed (red) line is the result 
for z = 0. One finds that TO3 has a minimal value around 0.015 eV, i.e. 7713 cannot vanish, and for 
z = it also has a maximal one around 0.027 eV. These two bounds can be found as well by using 
Eq. (|2ip . The non- vanishing of TO3 ~ \m ee \ agrees with the findings in the literature that \m ee \ is 
required to be larger than 0.01 eV, if neutrinos follow an inverted hierarchy [24]. Figure [2] shows 
that the relation in Eq. (|2ip . which is fulfilled to a good accuracy for z = 0, gives a lower bound for 
z 7^ in the tan Qn-frij, plane and no further constraints on the solar mixing angle can be derived. 
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Note that we used the best-fit value of the atmospheric mass squared difference for the dashed (red) 
line in Figure [2j Finally, we plot the Major ana phases <f>x and <p2 in Figure [3] against the lightest 
neutrino mass m^. As one can see, the phase 4>\ (blue/darker gray) varies between ir/8 and 7tt/8, 
while <p2 (green/lighter gray) either lies in the interval [0, 7r/8] or [7tt/8,tt] for small values of 1713, 
i- e - m 3 0.06 eV. The dashed (red) lines indicate again the value of 4>\ and 4>2 achieved in the 
limit z = 0. As the difference <j>i — <p2 of the two Majorana phases is the only quantity which can 
be realistically determined by future experiments [22] through 

\m ee \ = \mi cos 2 6 12 e 2i + m 2 sin 2 12 \ , (29) 

we also plot (pi — <p2 against 7713 in Figure [H This plot shows that the phase difference has to lie in 
the rather narrow ranges [— 3-7r/4, — vr/2] or [tt/2, 37r/4] for small values of m^. As one can see, the 
deviations from \<j>\ — 02 1 = 7r/2 (z = case) become larger for larger values of 1713. 



3.3 Flavon Superpotential 

In the following we discuss the flavon superpotential and show that the VEV structure assumed 
in (Eq.([7]) and) Eq.([8]) naturally arises, as does the spontaneous CP violation. In constructing the 
superpotential we work along the lines of [12]. For this purpose, we generalize -R-parity to a U(1)r 
symmetry under which the "matter fields" transform with charge +1, the fields h u and hd and 
the flavons are uncharged and another type of fields, the driving fields, have charge +2. These 
fields transform trivially under the Standard Model gauge group, but non-trivially under the flavor 
symmetry. The set needed for constructing the potential consists of Xe ~ (!l> a;4 )) a ° ~ (14)^) 
and Xu ~ (il> w ) under {D^Z§). Since all terms of the superpotential have to have U(1)r charge 
+2, the driving fields cannot couple to the fermions and can only appear linearly in the flavon 
superpotential. The renormalizable D4 x Z5 invariant superpotential for flavons and driving fields 
reads 

w f = ax e xl + bx° e <pl (30) 
+ c a {il -iD+dxlii^ + exlrf + f xl xl ■ 

Assuming that the flavons acquire their VEVs in the supersymmetric limit we can use the F-terms 
of the driving fields to determine the vacuum structure of the flavons. The equations 

a>xl + b<pl = Q , (31a) 

c(^-Vl)=0, (31b) 

d^^ + etpl + fxl = 0, (31c) 

result in 

fc) = ± i\/ffe), = M = ±! ^<M<M±Z<^Z (32 ) 

which can be re- written as 

, . fb .. . , , . I dj^j^ + e^ 

w e = ±1 y ~U e , {lpi) = ±V, W = ±lJ . (33) 



uw 



f 



dX°e 

dwj 

dwf 
~dxl 
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Note that the VEVs (ip e ) = u e , {1P2) = v and {ip v ) = u are unconstrained by the potential. Note 
further that the choice of sign in all cases is independent in Eq. (|32h and Eq. (|33p . For the discussion 
of the preserved subgroup structure it is anyway only relevant whether (tpi) = (^2) or (^1) = — (^2)- 
For (^i) = (^2) as used in Eq.([8]) we conserve a subgroup Z2 of D4 generated by B, whereas the 
relation (^i) = — {1^2) indicates that the Z2 subgroup generated by BA 2 is left unbroken. This Z2 
group is also not a subgroup of the D2 group conserved in the charged lepton sector. Thus, the 
subgroups of the charged lepton and the neutrino sector will be misaligned in both cases. In this 
paper we only consider the case of (ipi) = (^2) = v. Eq. (j33j) shows then that the VEVs w e and w 
necessarily have to be imaginary, so that CP is spontaneously violated, if the parameters a,...,f 
and the VEVs u e , v and u are chosen as positive. 

We remark that due to the U(1)r symmetry a /i-term fih u hd is forbidden in our model and 
has to be generated by some other mechanism. This feature is shared by all models using a U(1)r 
symmetry. In the derivation of Eq. (|3ip terms of the form Xu^uhd which couple a driving field 
to the MSSM Higgs fields can be safely neglected. They also cannot induce a //-term, since only 
vanishing VEVs are allowed for the driving fields, if the parameters a, ...,/ and the flavon VEVs 
are non-zero, as it is in our case. Finally, note that we find flat directions in this potential in the 
case of spontaneous CP violation under discussion here. These are however expected to be lifted 
by the inclusion of the NLO corrections, see Section 14.2^ as well as through soft supersymmetry 
breaking terms. Such nearly flat directions might be of interest for inflation [25]. 



4 Next-to-Leading Order Corrections 

In order to determine how our results are corrected at NLO, we take into account the effects of 
operators which are suppressed by one more power of the cutoff scale A compared to the LO. Such 
contributions to the fermion masses include two instead of only one flavon. In the flavon superpo- 
tential we add terms consisting of one driving field and three flavons. It turns out that there are 
actually no contributions to the fermion masses from two-flavon insertions due to the Z5 symmetry. 
Hence, the only NLO corrections we need to consider are those of the flavon superpotential, which 
lead to a shift in the flavon VEVs parameterized as 

(Xe) = w e + 5w e , (xu)=w + 5w and (ipi) = v + 5v . (34) 

The VEVs {(p e ) = u e , (ip u ) = u and (V'2) = v which are not determined at LO remain unconstrained 
also at NLO. The natural size of the VEV shifts is 

{VEV A 2 . (35) 



VEV 



As will be discussed in Section T4.21 the shifts 5w and 5w e are in general complex, whereas the shift 
5v in the VEV (ipi) is real for this type of spontaneous CP violation. 



4.1 Fermion Masses 

The VEV shifts induce corrections to the lepton mass matrices given in Eq.([9]) when the shifted 
VEVs are inserted into the LO terms, see Eq.@. In case of the charged lepton masses only the 
VEV of Xe is shifted. Such a shift is however not relevant, since it can be absorbed into the Yukawa 
couplings yf and y|- Especially, Ui is still given by Eq. (|10|) . The form of the neutrino mass 

6 These then become complex which however does not affect our results. 
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matrix is changed through the shifts of the VEVs into 

„,2 / yi(w + 5w) y 2 v 



M = — 
A 2 



y 2 (?; + £«) 
2/2^ y 3 w 2/4(10 + 5w) 

y 2 (v + 5v) yi(w + 5w) y 3 u 



(36) 



Note that 5w cannot be simply absorbed into uu, since 8w is complex, whereas w is imaginary. In 
the charged lepton mass basis the matrix in Eq. (|36p reads 



(37) 



yi (w + 5w) y 2 (v + e i7T ^5v / V5) y 2 {v + e - i7r / 4 5t; /y/5) 

K = j± [ y 2 (v + e^Sv/VS) y 4 (^ + <H 

y 2 (i> + e - " 1 "/ 4 ^/^) y^u y A (w + 5w) 



To evaluate the shifts in the neutrino masses and to discuss the deviations of the mixing angles 
from their LO values, especially #13 from zero and #23 from maximal, we parameterize the Majorana 
neutrino mass matrix as 



2 / is(l + ae) t(l + e in / 4 e) t (1 + e"^/ 4 e) 
M' u = ^jl t(l + e W4 e ) i x (l + ae) z 

A 1 t(l + e-^/ 4 e) z ix(l + ae) 



(38) 



with s, t, x and z as given in Eq. (|17|) and0 



$ w ■ 1 1 $ v 2 

q e = — , a = a r + 1 ^ ana e = — = — « A U.U4 



(39) 



The neutrino masses and mixing parameters resulting from Eq. (|38p can then be calculated in an 
expansion in the small parameter e. We observe that the mass shift of m 2 would vanish for Sw 
being zero. Its explicit form is 



(40) 



with (m^°) 2 given in Eq. (|18p . Similarly, the masses m 2 and m| undergo shifts proportional to e. 
A simple expression can however only be found for the sum m\ + m 2 



K LO ) 2 + (mf ) 2 = K°) 2 + {m\°f + 2 (|) 2 g) a (2 V2t 2 + a r ( S 2 + x 2 ) -a i xz)e (41) 

(m^) 2 can be found in Eq, (|18p . The mixing angle #13 no longer vanishes and we find 

t x 



sin 013 



For 6*23 we get 



t 2 + (s — x)x — z 2 
tan0 2 3« 1 + V2 



e . 



x z 



t 2 + {s — x)x — z 2 

The deviation from maximal atmospheric mixing can also be expressed through 



cos26> 23 | ~ V2 



X z 



t 2 + (s — x)x — z 2 



V2 



sm 013 



(42) 
(43) 

(44) 



7 We assume that e is positive. 
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From both formulae one can deduce that in the case z = the corrections to maximal atmospheric 
mixing are not of the order e, but only arise at 0(e 2 ). Contrary to this #13 still receives corrections 
of order e, if z = 0. The solar mixing angle 6*12 which is not fixed to a precise value in this model also 
gets corrections of order e. We note that the smallness of \s — x\ and z, required by the smallness 
of Am^i, might lead to a disturbance of the expansion in the parameter e. 

A correlation between cos 26*23 an d sin #13 depending only on physical quantities, Am|, and 
not on the parameters of the model, s, t, can be obtained by an analytic consideration which is 
done analogously to the study performed in [7]. Clearly, the matrix in Eq, (|37p is no longer fj, — r 
symmetric, however we find the following remnants of this symmetry 

(MX = (^X and (M' V )^ = (M' V ) TT . (45) 

Eq. (|45p shows that fJ, — T symmetry is only broken by phases, but not by the absolute values of the 
matrix elements. This leads to 

= |(MXl 2 + l«WI 2 - \(K)er\ 2 ~ \(M' U ) TT \ 2 (46) 

3 

= {KK^ ~ (KM'J)rr = J2 m2 M U MNs)„j\ 2 ~ \{U M N S )r 3 \ 2 ) 

i=i 

= ^(sin 2 612 — sin 2 #13 cos 2 6*12) m\ + (cos 2 #12 — sin 2 #13 sin 2 #12) m 2 — cos 2 #13 rri^j 005(26*23) 
— Am^i sin 26*12 sin 26*23 cos § sin $13 • (47) 

Since sin 6*13 ~ 0(e) and cos 26*23 ~ C( e ) is already known, we can linearize this equation and 
obtain (using best-fit values for the physical quantities and the fact that neutrinos have an inverted 
hierarchy in this model) 



Am|i sin 26*x2 
Am| + Am 2 .! sin 2 6*12 



cos 26*23 ~ 5 — — n — r~o cos 5 sin #13 0.03 cos 5 sin 6*13 . (48) 



Eq. (|48p can be used to estimate the largest possible deviation from maximal mixing. For sin 6*13 
being at its 2 a limit of 0.2 and | cos<5| = 1, | cos 2^23 1 still has to be less than 6 x 10~ 3 which is 
well within the la error. Finally, we note that Eq. (|48p must be consistent with Eq. (|44[) and thus 
we again find that z ought to be small. 

4.2 Flavon Superpotential 

The corrections to the flavon superpotential stem from terms involving one driving field and three 
flavons. These terms are non-renormalizable and suppressed by the cutoff scale A. We find 

An// = I X° e Xl + k f Xl X, Vl + § X° e X»M2 + ^X° eV >» (V? + V> 2 2 ) (49) 



Assuming that CP is only spontaneously violated forces all fcj to be real. We calculate the F-terms 
of Wf + Awf for the driving fields using that the VEVs can be parameterized as 

(Xe) = w e + 5w e , (xu) = w + 5w and (^i) = v + 5v . (50) 
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The VEVs (<p e ) = u e , ((p v ) = u and (^2) = v are not determined at LO. We assume that only 
terms containing up to one VEV shift or the suppression factor 1/A, but not both, are relevant. 
The F-terms then lead to 

2 a w e 5w e + — (ki w 3 + u 2 w + k% v 2 w + 2 k^ uv 2 ) = , (51a) 
2cvSv + ^j(k 5 w 2 + k 6 ul) = 0, (51b) 
dv5v + 2fw5w + ^(k 7 wl + k 8 u 2 e ) =0 . (51c) 

Here we have chosen the solutions with + in Eq. (|33p . The explicit form of the shifts reads 

1 u 

5v = -- e - (k 5 w\ + k 6 u 2 e ) , (52a) 

2c v A 
I 1 

5w = - — - — - (d (k 5 wl + k^u 2 ) u e - 2 c (kr w 2 + kg u 2 e ) w e ) , (52b) 
4c/ w A 

5w e = r (k\ w 3 + k2 u 2 w + k% v 2 w + 2 k^u v 2 ) . (52c) 

2 a w e A 

As one can see, for our type of spontaneous CP violation 8v is real, whereas 5w e and 5w turn out 
to be complex in general. As can be read off from Eq. (|52p all shifts are generically of order 

^^~A 2 for VEV ~ A 2 A. (53) 
VEV 

Finally, note that the free parameters (<p e ) = u e , {tp u ) = u and (^2) = v are still undetermined. 



5 Summary and Outlook 

We constructed a supersymmeterized version of the D4 model by Grimus and Lavoura [1]. For this 
purpose, we replaced the Higgs doublets transforming under the flavor group D4 by gauge singlets. 
We also enlarged the auxiliary symmetry which separates the different flavor breaking sectors from 
2,(aux) ^ 2^ The simplest supersymmeterized -D4 model does not contain right-handed neutrinos, 
but neutrinos get masses through the operator lh u lh u /A. Apart from these slight changes the 
model is essentially the same as the one by GL, since we also generate maximal atmospheric 
mixing and vanishing #13 through the fact that -D4 is broken to D2 in the charged lepton and to 
Z2 in the neutrino sector. The crucial issue of the vacuum alignment is elegantly solved through 
an appropriately constructed flavon potential. We performed a phenomenological analysis under 
the assumption of a certain type of spontaneous CP violation suggested by the minimization of 
the potential. As a result, the neutrinos have to have an inverted hierarchy. The quantity \m ee \, 
measured in Ov/3/3 decay, is almost equal to the lightest neutrino mass m^. Furthermore, we found 
that 771,3 cannot vanish and has a lower bound around 0.015eV. The Majorana phases 4>i and (f>2 are 
restricted to a certain range at least for small 771,3 • 111 contrast to that the solar mixing angle #12 can 
take all values allowed by experiments. We also analyzed the NLO terms in this model and showed 
that they only induce shifts in the VEVs of the flavons, but no additional terms in the Yukawa 
sector appear. The shifts yield deviations from the LO results, 613 = and 623 = it/ 4. Comparing 
these deviations we see that although both of them could in principle be of order e ~ A 2 ~ 0.04, 
the smallness of the parameter z, necessary to arrive at mass squared differences and 612 within 
the 2a ranges, leads to the fact that 623 is much closer to 7r/4 than #13 to zero. 
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The supersymmeterization of the D4 model has to be regarded as a step towards a grand unified 
model with a dihedral flavor symmetry for two reasons: (a) low scale supersymmetry elegantly cures 
the hierarchy problem and easily allows the gauge couplings to be unified at 10 16 GeV and (b) the 
replacement of the Higgs doublets transforming under the flavor group by flavons is important for 
disentangling the breaking of flavor and gauge groups. However, another essential feature of a 
unified theory is that the lepton sector cannot be discussed without also considering quarks. So, 
one of the major challenges to tackle is the question how to implement the quark masses and their 
mixings in a model with a dihedral symmetry. In the recent past models have been presented which 
are able to predict the Cabibbo angle with the help of the flavor group Dj and -D14 [8-10]. Thus, it 
is interesting to search for a way to combine these models and to find a (probably larger) dihedral 
symmetry leading to the same results, which we get from a D4 flavor group in the lepton sector 
and from a D7 or D14 group in the quark sector. 
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A Connection to the Group Basis chosen in [1] 

Note that replacing the left-handed fields I by Ufl, with U[ given in Eq.ljlOp. is equivalent to 
changing the basis in which the generators A and B are given for the two-dimensional representation 
of -D4. Since also the second and third generation of the right-handed charged leptons form a doublet 
under D4, we also have to transform them to show that this corresponds to a change of the generator 
basis of the D4 doublet. By calculating the matrix U e c which diagonalizes M} Mi one finds that U e c 
equals Uf P where P is a diagonal phase matrix. This matrix P induces an unphysical rephasing 
of the right-handed fields to keep Jjj Mi U e c a diagonal matrix with positive entries. The change of 
basis (induced by the unitary matrix U*) leads to real generators A and B of the form 

This coincides with the basis chosen in [1], if A is identified with the product hg and B with the 

generator h. The identifications for the singlets are the following: |- corresponds to l_i 

to I4, 1 |_ to I3 and 1 to 1_2- In [1] the charged lepton mass matrix is diagonal without any 

further transformation and #13 = and #23 maximal can be directly read off from the neutrino 
mass matrix. This is the same in our case, if we go into the primed basis, see M' v in Eq.p3]). 

B Importance of Mismatch of Subgroups 

To elucidate the reason why the two subgroups preserved in the charged lepton and the neutrino 
sector have to be different, i.e. the subgroup present in the neutrino sector should not be a 
subgroup of the D2 group of the charged lepton sector, observe that M\ m} as well as M v Ml for 
Mi and M v given in Eq.® can be written in the following form 

MiM] = Bie~ il3i d Die**** i = l,v. (55) 

\ Bie -< Ae-^ j Ci ) 
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This form is achieved for M; (M u ) as long as at least a Z 2 group, originating from B A m , is conserved 
in the charged lepton (neutrino) sector. A matrix of this type is diagonalized through 

jPi n 


e" 



Ui 



-i <Pi) 



U max U l2 (0i)U(ai) 



(56) 



where 








( cos 9i 


sin 9{ 





u 12 (e 


H 


— sin 6i 


cos 6i 









V o 





1 



U(al) 



V 












(57) 

In [9, 10] it has been shown that the quantities <f>{ and j are related to the group theoretical indices 
of the flavor symmetry, j is the representation index of the doublet under which two of the three 
left-handed lepton generations transform. Thus, it is the same for charged leptons and neutrinos. 
4>i can be expressed as 

(58) 



rrii 



n 



where n is the index of the group D n and mi r v \ the index of the preserved subgroup in the charged 
lepton (neutrino) sector that has a generator of the form B A" 1 ' M . mi ( v \ is an integer number 
between zero and n — 1. The parameters Ai, Di and the phase are real functions of the 
matrix entries of MiMl, whose proper form is not needed here. The phases at are irrelevant for 
the diagonalization of Mj M.\ , but are necessary for the diagonalization of the neutrino mass matrix 
M v alone. The angle 9% can be expressed through the parameters Ai, Di as follows 

2^2 Bi 



tan2#j = 

The general form of the MNS matrix is then 



Ci + Di — Ai 



(59) 



Umns = Uf U* 



U(a{)U? 2 ( 



T 




U 12 (9 U )U(- 



1 

e 

(60) 

This form already shows that it is essential to have a non-trivial phase e ^ v '' } in order to 
guarantee that the maximal mixing in the 2 — 3 sector is not cancelled. For the third column of 
Umns, which determines the mixing angles #13 and #23 > we find 



\{U M Ns)ez\ = I sin((0; -<f> v )]/2) sin 0/| , \{U M ns)^\ = I sin((^ -(f> v )j/2) cos6»;| 
\(Umns)t3\ = I cos((0/ - <f> v ) j/2)| . 



(61) 



Using that we preserve a Z 2 symmetry generated by B in the neutrino sector and a D 2 group 
generated by BA (according to our convention for the generators of the group D 2 introduced in 
Section [2]) in the charged lepton sector, gives for <\> v and 

0„=O and <k = \- (62) 

j is trivially one, since D4 only contains one irreducible two-dimensional representation 2. As the 
elements (1,/s) and (k, 1) with k = 2, 3 in Mi vanish, see Eq.Q, the parameter Bi in Eq. (|55p is 
zero (and also 0i = 0) and thus 9i = as well according to Eq, (|59p . This results in 

\{UMNS)e3,\ = , \{UmNs)/i3\ = \(UmNs)t3\ = (63) 
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giving maximal atmospheric mixing and vanishing #13. A few things are interesting to notice: 
In principle four different cases might occur. These arise from whether the subgroups D2 and 
Z2 contain the same element BA m or not and from whether the D2 subgroup is unbroken in the 
charged lepton sector or only a Z2 subgroup is preserved. The first issue determines whether mi 
equals m v or not, i.e. whether \4>i — <f> v \ is zero or not. The second one is responsible for (non-)zero 
0\. We can see from Eq. (161j) that for no mismatch of the subgroups #13 as well as #23 vanish, in 
contrast to what is observed in nature. So the mismatch of the two subgroups is necessary. If 
61 is zero, i.e. the subgroup present in the charged lepton sector is a D2 group, #13 = and 621, 
maximal follow. If however only a smaller Z2 group is present in the charged lepton sector, neither 
013 being zero nor #23 being maximal holds. Then only the MNS matrix element \{Umns)t?\ is 
fixed by group theory. 

Finally, the matrix Ui given in Eq. ljlOp equals the matrix shown in Eq. (|56p . if we additionally 
set the phases to a[ = 0, a l 2 = vr/4 and a 3 = 37r/4. 

One might ask the question what actually determines the size of the solar mixing angle 612 in 
this context. For Q\ = we find from Eq. (|60p that 

\(UMNs)ei\ = \cos9 u \ and | (U M Ns)e2\ = I sin V \ (64) 

which shows that 612 is given by 9 U . Since this angle would vanish, if a D2 group instead of a Z2 
group (with generator BA m ) was present in the neutrino sector, one might interpret the size of the 
solar mixing angle as hint to how strongly a D2 group is broken in the neutrino sector. 
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